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Abstract 

The paper is concerned with a system consisting of two coupled nonlinear parabolic equations with a cross-diffusion 
term, where the solutions at positive times define the initial states. The equations arise as steady state equations of 
an age-structured predator-prey system with spatial dispersion. Based on unilateral global bifurcation methods for 
Fredholm operators and on maximal regularity for parabolic equations, global bifurcation of positive solutions is 
derived. 
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1. Introduction 

Consider the situation that an age-structured prey and an age-structured predator population inhabit the same 
spatial region 0,, that the individuals of both populations undergo spatial fluctuation and that the predator population 
exerts a repulsive pressure on the prey population. If u = u(t, a, x) > and v = v(t, a,x) > denote the density 
of the prey and the predator population, respectively, at time t > 0, age a £ [0, a m ) for a maximal age a m > 0, and 
spatial position x £ 17, a simple model reads 

dtu + d a u — A x ((Si + jv) u) — —aiu 2 — a2uv , t>0, a £ (0, a m ) , x £ Q, , (IT) 
d t v + d a v - 6 2 A x v = -/3 1 v 2 +p 2 vu, t>0, oe(0,a m ), x e Q . (1.2) 

These equations are subject to the nonlocal age boundary conditions 

u(t,0,x)= / B 1 (a)u(t,a, x) da , i > , i6fi, (1.3) 
Jo 

v(t,0,x)= / B 2 (a)v(t,a,x) da , t>0, 160, (1.4) 
Jo 

the spatial boundary conditions 

u(t, a, x)=0, t>0, a £ (0, a m ) , x £ dSl , (1.5) 
v(t,a,x)=0, t>0, a£(0,a m ), x £ 80. , (1.6) 

and are supplemented with time initial conditions. The A x -term in ( ll.lt describes spatial movement of prey individ- 
uals. Besides intrinsic dispersion with coefficient 5% > 0, it reflects an increase of the dispersive force on the prey by 
repulsive interference with an increase of the predator population. Here, 7 > is the predator population pressure 
coefficient. We refer to 112411 for a derivation of such kind of models (without age-structure). The right hand sides of 
(ll.lt and (11.21 ) take into account intra- and inter-specific interactions of the two populations with positive coefficients 
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a\, 0.2,01,132 > 0. Equations ( 11.5b and ( 11.61 ) describe creation of new individuals with nonnegative birth rates B\ 
and i?2- The reader is referred to [32] and the references therein for further information about linear and nonlinear 
age-structured population equations with (and without) spatial dispersal. Although such models have a long history, 
there does not seem to be much literature about equations with nonlinear diffusion. The well-posedness of equations 
( II . lb - dT~6b might be derived within the general framework presented in 112.811 though the results therein do not directly 
apply due to the cross-diffusion term. 

We shall remark that the model considered herein is a rather simple biological model and there may be more 
accurate models, e.g. with different maximal ages, nonlocal dependences in the reaction terms etc. The aim of the 
present paper is thus rather to provide a mathematical framework to treat parabolic equations with cross-diffusion and 
nonlocal initial conditions. Namely, the present paper is dedicated to positive equilibrium (i.e. time-independent) 
solutions to ( ll.ll )-( fL6l l which shall be established based on global bifurcation methods. We write the birth rates in 
the form B\(a) — ?/6i(a) and B2(a) = £62(0), where 61, 62 are some fixed birth profiles and the parameters rj, 
£ measuring the intensities of the fertility shall serve as bifurcation parameters. Aiming at a simple and compact 
notation, we let 6 := bi = 62 and 61 = 82 = 1. We emphasize that these simplifications are made merely for the 
sake of readability and do not impact in any way on the mathematical analysis to follow. To shorten notation further, 
we shall agree upon the following convention: given a function defined on J := [0, a m ] and denoted by a lower case 
letter, say u or v, we use the corresponding capital letter U or V to denote its age-integral with weight b, i.e., 

U := / b(a)u(a)da, V := / b(a)v(a)da. (1.7) 
Jo Jo 

Depending on the values of the parameters 7/ and £, we are looking for nonnegative and nontrivial functions u — 
u(a, x) and v — v(a, x) satisfying the nonlinear parabolic system with cross-diffusion 

d a u — A x ((l + jv)uj — — aiu 2 — 0L2UV , a S (0, a m ) , x e f2 , (1-8) 
d a v - A x v = -Pxv 2 + f3 2 vu , ae(0,a TO ), x e fl , (1.9) 



subject to the nonlocal initial conditions 



and spatial Dirichlet boundary conditions 



u(0,x)=r)U, xen, (1.10) 
1,(0,2;)=^, (1.11) 



u(a,x)—0, dg(0,a m ), x £ <9f2 , (1-12) 
v(a,x)=0, ae(0,a m ), x S dil . (1.13) 

Clearly, of particular interest are coexistence states, that is, solutions (u, v) with both components u and v nontrivial 
and nonnegative. 

Based on bifurcation techniques, existence of equilibrium solutions was established in 1 3(J 31 ] for similar systems 



with linear diffusion and in [26, 27, 29] for a single equation with nonlinear diffusion. Except for [29], the bifurcation 
parameter was also chosen to be a measure for the intensity of the fertility as above. Prior to the just cited papers, 
bifurcation methods were used in [9] to derive positive equilibrium solutions for a single equations with linear dif- 



fusion. We also refer to [ 15] where the large time behavior of population dynamics with age-dependence and linear 
spatial diffusion was analyzed. 

More attention attracted than the age-structured parabolic equations so far have related elliptic systems of the form 

— A x [(l + 777)?/] = u(fi — u ± cv) , — A x v = v(X — v ± bu) 

in fi subject to Dirichlet conditions on dfl with positive constants /1, A, 6, c; both for the case of linear diffusion 
(e.g., s ee [Hi 0, S OS [l8ll and the references therein) and the case with cross-diffusion 7 > (e.g., see 



10, 1 U, [131, HJ, [23J and the references therein). It is worthwhile to remark that for such elliptic equations with cross- 



diffusion, the transformation z := (1 + 777)?/ leads to a semilinear elliptic system for z and v, for which, when written 
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in the form (z, v) — S(z, v), the corresponding solution operator S enjoys suitable compactness properties. More- 
over, for such a system, positivity and a-priori bounds of solutions may be derived using the maximum principle, e.g. 
ifiol [l3 . 231. For the nonlocal parabolic equations dl.7t -( [TT3T l under consideration herein, however, a correspond- 
ing transformation still yields nonlinear second order terms with respect to the spatial variable (or first order time 
derivatives). Consequently, the underlying solution operator does not enjoy similar compactness properties. As we 
are interested in global bifurcation of positive solutions, it is thus not clear how to apply directly the unilateral global 
bifurcation methods of |T^, 21] as in \ 3(&, which require compactness of the underlying solution operators. To over- 
come this deficiency, we shall invoke recent results of Shi & Wang 112311 on unilateral global bifurcation for Fredholm 
operators, which are based on Lopez-Gomez's interpretation 11611 of the global alternative of Rabinowitz 112 111 and on 
the global bifurcation results of Pejsachowicz & Rabier lfl9ll . These results yield a (global) continuum of positive 
solutions. We shall also point out that, due to the cross-diffusion term, no comparison principle in the spirit of l30l 
Lem.3.2] is available. In some cases, the lack of such a comparison principle prevents us from determining which of 
the possible alternatives the constructed continuum of coexistence states satisfies (see Theorem 12. 3 1 and Theorem l2.4l 
below). 



2. Main Results 

We suppose throughout this paper that the birth profile 

b G ((0, a m j) with b(a) > for a near a m (2.1) 

is normalized such that 

b(a)e- Xia da = 1 , (2.2) 

where Ai > denotes the principal eigenvalue of — on £1 subject to Dirichlet boundary conditions on cMl. Before 
stating our results on coexistence states in more detail, we first recall some auxiliary results from 13011 about semi- 
trivial states, that is, about solutions (it, v) with one vanishing component. Taking for instance v = in dl.8t - dl.9t 
and using convention dl.7t . we obtain the reduced problem 

d a u — A x u = —a.\u 2 , u(0, •) = rjU , (2.3) 

subject to Dirichlet boundary conditions. Introducing, for q G (n + 2, oo) fixed, the solution space 

W q := L q ((0,a m ),Wl D (Q)) nW^((0,a m ),L g (n)) > 

where 

Wl D {VL) := {u G W^{VL) ■ u = OondO} , 
and letting W+ denote its positive cone, the following result was proved in B30l Thm.2.1, Cor.3.3, Lem.3.6, Cor.3.7]: 

Theorem 2.1. Suppose d2.lt and ( 12.2b . For each rj > 1 there is a unique solution u v G \ {0} to equation 
d2.3t . The mapping (rj i— > u n ) belongs to C°°((l, oo), W q ) and satisfies \\un\\w t as V ~ * ^ and \\u 1l \\^ q — > oo 
as rj —¥ oo. For rj > 1 and a G J, the derivative d rt u rl (a) is a (strictly) positive function on f2. If r\\ > r\i, then 
u Vl > Ujj 2 - Finally, if r/ < 1, then d2.3t has no solution in \ {0}. 

Given £ > 1, we let G W+ \ {0} denote the unique solution to the corresponding equation for v when taking 
u = in ( 1 1 .8b - dTT9b . Though we shall work in the solution space W q , we remark that solutions to dl.7l )-( fTTT3b . so in 
particular u v and v%, are smooth with respect to both a € J and x G fi (see Lemma [3~3l below). 

2.1. Bifurcation with Respect to the Parameter rj 

We first shall keep £ fixed and regard rj as bifurcation parameter. We thus write (77, u, v) for solutions to (11.7b - 
( 11.1 3b with u, v belonging to W+. Theorem l2.1l entails. for any £ > 0, the semi-trivial branch 

Bi := {(77, u v , 0) ; 7/ > 1} C K+ x (W+ \ {0}) x W+ (2.4) 
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of solutions. For £ > 1, an additional semi-trivial branch 

03 2 := {(??, 0, v 6 ) ■, v > 0} C M+ x W+ x (W+ \ {0}) (2.5) 
exists. In the latter case, a continuum of coexistence states bifurcates from 032: 

Theorem 2.2. Suppose d2.1t -( l2~2l and let £ > 1. There exists 770 := ?7o(£) > smc/i f/?af an unbounded continuum £ 
emanates from (?7o, 0, v%) £ 032, a«<^ 

£ \ {(770, 0, c M+ x (W+ \ {0}) x (W+ \ {0}) 

consists of coexistence states (77, it, v) to dl.7t - dl.13t . Afear f/ie bifurcation point (770, 0, Wg), £ ;s a continuous curve. 
There is no other bifurcation point on 032 or on 03i to positive coexistence states. 

The value of 770 (£) corresponding to the bifurcation point is determined in ( I3.30t . 

We turn to the case £ < 1, which is more involved. Recall that 03 1 is the only semi-trivial branch of solutions in 
this case. 

Theorem 2.3. Suppose d2.lt - d2.2t . There is 5 G [0, 1) with the property that, given £ G (5, 1), there exists 
Vi := T /i (0 > 1 swc/i f/ia? a continuum & emanates from (771 , u,^ ,0) G 03i, ant/ 

6 \ {(77x, Ut?1 , 0)} c R+ x (W+ \ {0}) x (W+ \ {0}) 

consists of coexistence states (77, u, v) of dl.7t - dl.13t . 77ze continuum & 

( i) is unbounded, or 

(ii) connects (r]i,u rn ,0) to a solution (77*. n. v) of dl.7t - dl.13t with n* G (0, 1) and u, v G \ {0}. 

Near the bifurcation point (r)i , u m , 0), © is a continuous curve. There is no other bifurcation point on 03 1 to positive 
coexistence states. 

The values of 5 and 7/1 (£) are given in ( 14. 5t and ( 14. 6t . respectively. Note that alternative (i) above always occurs 
if cross-diffusion is not taken into account, i.e. if 7 = 0, see jjll Thm.1.5]. Alternative (ii) stems from a technical 
condition, and we conjecture that alternative (i) is the generic case also for the present situation where 7 > 0. However, 
as occurrence of alternative (ii) is not ruled out by our analysis, Theorem l2.3l is rather a local bifurcation result in this 
regard. 

2.2. Bifurcation with Respect to the Parameter £ 

Since the cross-diffusion term involves the predator density v, bifurcation from semi-trivial solution branches is 
more intricate when regarding £ (the measure of the predator fertility) as parameter and keeping 77 fixed. We shall 
only consider the situation 77 > 1, the reason being explained in Remark IBTTI We now write (£, u, v) for solutions to 
dl.7t - dl.13t . Then, there are two semi-trivial branches of solutions 

Ti := {&0,tO; £ > 1} , T 2 := 0) ; £ > 0} . 

Also in this case, bifurcation from the semi-trivial branch T2 occurs: 

Theorem 2.4. Suppose d2.lt - d2.2t and let 77 > 1. There exists £0 : — £o(??) £ (0, 1) such that a continuum *K of 
solutions to dl.7t -( [TTT3l emanates from (£0, U^, 0) G T2 satisfying the alternatives 

(i) m \ {(£0, u v , 0)} is unboundedin M+ x (W+ \ {0}) x (W+ \ {0}), or 

(ii) connects T2 with %\. 

Except for the bifurcation point(s), ?t consists exclusively of coexistence states. Near the bifurcation point (£0, u v , 0), 
9\ is a continuous curve. There is no other bifurcation point on %2 or on Ti to positive coexistence states. 
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The values of (,o(i]) and £1 = £1(77) > 1, corresponding to the connection point (£1, 0, v^) £ Ti if alternative (ii) 
occurs, are determined in ( 15. U and ( 15.21 ). respectively. Given the assumptions of this theorem, if 7 = 0, one can show 
that alternative (ii) occurs if 77 is less than a certain value or if additional assumptions on the birth rates and on otj, f3j 
are imposed, see IDOL Thm.2.2]. 

The outline of the paper is as follows: The next section is dedicated to the proof of Theorem l2.2l In Subsection l3.11 
we first introduce some notation, and in Subsection l3.21 we provide the necessary auxiliary results needed to perform 
the actual proof of Theorem l2.2l in Subsection [33] Sections [4] and [5] are dedicated to the proofs of Theorem 12. 3 1 and 
Theorem |231 respectively. As these proofs are along the lines of the proof of Theorem 12. 21 these sections are kept 
rather short. 

3. Proof of Theoreml£2l 

3.1. Notations and Preliminaries 

Given two Banach spaces F\ and Fq, we shall use the notation C(Fi, Fq) for the set of bounded linear operators 
and IC(Fi,Fq) for the set of compact linear operators between F\ and Fq. We set C(Fq) := C(Fq, Fq) and similarly 
K{Fq) := JC(Fq, Fq). The set of toplinear isomorphisms between F\ and Fq is denoted by Cis(Fi , Fq). 

Let —Ad denote the negative Laplacian on L q :— L q (Vt) subject to Dirichlet boundary conditions, that is, 
— An := —A x with domain W qD , where 

W£ D := W£ D (n) := {u e ; u = on dtl} 

for k > l/<7 and W£ D := W q (f2) for < k < l/q. As q > n + 2, we have W 2 ~ D 2/q -4 C^O) by the Sobolev 

embedding theorem. In particular, int(W qD 2 ' + ) =/= 0, that is, the interior of the positive cone of W^ D 2/q is not 
empty. Let 70 denote the trace operator defined by 70M := u(0) for u <G W 9 , which is well-defined owing to the 
embedding III.Thm.4.10.2] 

W q ^C([0,a m ],W 2 - D 2/q ) . 
In fact, we have, due to the interpolation inequality [2, I.Thm.2.1 1.1], 

W q ^C 1 - 1 '"-\[Q,a m ],W 2 q %) , 0<<?<1-1/V (3.1) 

In the following, we let W^" := \ {0}. Recall that an operator A S C(W q D , L q ) is said to have maximal 
L q -regularity provided 

{8 a + A 7o) 6 Cis(W q ,L q x W 2 ~ D 2/q ) , 

where h q := L q (J,L q ). Recall that — Ad has maximal L 9 -regularity. If A : J — > C(W q i D , L q ) for some k > 0, 
we write (Au)(a) :— A{a)u{a) for a 6 J and u : J — > W qD . The following perturbation result is proved in 1I20I 
Cor.3.4], which is also valid for Revalued functions, i.e., if L q is replaced by L q (il, M. k ) etc.: 

Lemma 3.1. Let A £ C(J, C{W q D , L q )) be such that Ala) has maximal L q - regularity for each a G J and let 
B e L q (J, C(W 2 - D 2/q , L q )). Then (d a + A + B, 7o ) £ £is{W q , h q x W 2 q - 2/q ). 

Given q > and h 6 C e (J, C(f2)), we let II^j (a, a), < a < a < a m , denote the unique parabolic evolution 
operator coiTesponding to —Ad + h € C e (j,jC(Wq D ,L q )), that is, z(a) = U^(a,a)$>, a £ (<r,a m ), defines the 
unique strong solution to 

d a z — Adz + hz = , « 6 [c, o m ) , z(cr) = $ , 
for any given a £ (0, a m ) and $ £ L q (see yl II.Cor.4.4.1]). Note that the evolution operator is positive, i.e. 

U [h] (a, <j)$ £ L+ , < a < a < a m , $ £ L+ . 
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Since — Ad has maximal L g -regularity, it follows from Lemma [37X1 that 

(d a -A D + h, 7o ) e Cis(W q , h q x W^ D Vv ) (3.2) 
and, in particular, U [h] (-, 0)$ G W q for $ e W^ /q . We set 

:= / b(a)n [h] (a,0)da . 
Jo 

Then S K-(W^ ^ ) owing to standard regularizing effects of the parabolic evolution operator EI^] and the 

compact embedding Wg )i3 ^ • Moreover, as II^] (a, 0) for a 6 (0, a m ) is strongly positive on W^~^ q by 

10, Cor. 13. 6], the same holds true for Hnj due to ( 12. It . that is, 

%]$Gint(<^ + ), $£< D 2/!,+ \{0}. (3.3) 

The corresponding spectral radius r(-ffrw) can thus be characterized according to the Krein-Rutman theorem (jl 
Thm.3.2] (see iflLem.3.1]): 

Lemma 3.2. For h G C e (J, C(f2)) with o > 0, f/ze spectral radius r(H^) > is a simple eigenvalue of 
with a corresponding eigenfunction belonging to 'mt(W 2 ^ q ' + ). It is the only eigenvalue of with a positive 
eigenfunction. Moreover, ifh and g both belong to C e (J, C(Q)) with g > h but g ^ h, then r(H\ g -\) < r(HnA). 



In particular, the normalization (12.21 i implies 

r(H m ) = l (3.4) 

since any positive eigenfunction of -Ad is an eigenfunction of Hr^ as well. Moreover, writing the solution to (12.31 l 
in the form u = H[ ai u] (■, 0)u(0), Theorem l2. ll together with Lemma l3~2l imply 

Vr{H [aiUii] ) = £_r(H [/3lV(] ) = 1 , r},£>l, (3.5) 
since ^(0),« ? (0)G(< D 2/9 ' + )\{0}. 
3.2. Auxiliary Results 



The aim is to apply the global bifurcation results of II23L Thm.4.3,Thm.4.4] in order to establish Theorem l2.2l We 
first provide the necessary tools. 

Let £ > 1 be fixed and let G W+ denote the solution to ( ll.7t - dl.13t with u = provided by Theorem 12.11 
Throughout we use the convention ( 11. 7t . Notice that for each a G J, the operator A\ (a), given by 

Ai(o)«:=-div je ((l+7W € (o))V a! u), u^W 2 qD , 

has maximal ^-regularity due to its divergence form, the positivity of G W 9 , (13.lt . and e.g. [2, I. Cor. 1.3.2, 
III.Ex.4.7.3, III.Thm.4.10.7]. Moreover, Ax G C(J, £{Wl D ,L q )) by (l3Jl . Noticing also that 

Bi(a)u := div x (uX7 x v s (a)) + a 2 v^(a)u , a G J , u G , 

defines an operator £?i 6 L g (J, C(W q ^^ q , L q )), it follows from Lemma |3~T1 that 

Tx := (9 a + A 1 +B 1 , lQ y 1 e C(h q x W^ D Vq ,W q ) (3.6) 

is well-defined. Observe that 

Af.u := (Ai = -A£»((l+7«f» + a 2 V£tt , it G . 
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From ( 13.11 ) and ( 13.21 ) we obtain 

T 2 := (d a -A D + 2/3 1 ^, 7o r 1 G C(h q x W 2 q ^'\ W g ) . 

To derive a bifurcation from a point (770,0, uj) G 5$2 for a suitable 770 = f?o(£)> we write (77, u, v) = (77, u, + w) 
for a solution to problem dl.7t -( fTT3l >, which is then equivalent to 



d a u - A D ((1 + jv^)u) = A D (jwu) - otiu 2 ~ a 2 u(v^ + w) , 
d a w - A D w = -/?iuj 2 - 2j3iV£W + + w)u , 

subject to 

u(0,x) = rjU(x) , x efl , 
w(0,x) =£W(x) , xefl. 

The solutions (77, u, w) to (13 .7b - (l3 . 1 Ob . in turn, are the zeros of the map F 

F(r),u,w) : 

where 



a G (0,a m ) 
a G (0,a m ) 



x 



u — Ti (Afl(7wu) — aiu 2 — a^uw , 77?/) 



16O, (3.7) 
xen, (3.8) 

(3.9) 
(3.10) 

¥ 9 x W 9 , defined by 
(3.11) 



W q := {w G W 9 ; uj(o, x) > -1/27 , a G J , x G ft} 

is an open subset of W g owing to ( 13.11 ). Clearly, F is smooth and -F(t7, 0, 0) = for 7/ G 
derivatives at (77, u, w) we compute 



As for the Frechet 



F,, 



){r],u,w)[(t>,ij)\ 



cf) — Ti (Ad i'fipu) + Ac(7t(j0) — 2a\ucf) ~ a2W(p — ct2Uip , 77$) 
i>-T 2 (-2Piwtl) + Mu + l32(ve + w)<l>, £#) 



and 



-^ri,(u,w) (.V: 



-T x (0, $) 




(3.12) 



(3.13) 



for (0, t/j) G W 9 x W g . The choice of W 9 as solution space is to have a suitable functional setting to work with in 
the framework of maximal regularity. However, as it is needed later on, we note that solutions to ( ll.7t - dl.13t . i.e. to 
(I3.7t -( l3.10t . are smooth. The proof is a bootstrapping argument which we provide for the reader's ease. 

Lemma 3.3. If (r]j, Uj, vj) is a bounded sequence in K x W g x W g of solutions to (II .7b - (fTTT~3~b . then (uf) and (i)j) 
are bounded in C e (J, C 2+£ (ft)) D C 1+£ (J, C e (tt)) for some e > 0. 

Proof. To stick with the notation of |2|, let (E, A) :— (L q , —Ad) and let [(E a , A a ); a > 0] be the corresponding 
interpolation scale induced by the real interpolation functors (•, -) a .q- Putting 



F := E 1 _ 1 / q = D q , Fx \= 

it follows from [2, V.Thm.2.1.3] that the Fq -realization of Ad, again denoted by Ad, has domain F\ and is the 
generator of an analytic semigroup {e aA ° ; a > 0} on Fq. Thus, 

\\e aAD \\ c{F ^ Fv) <c a^ , aGJ\{0}, G (0, 1) , (3.14) 

where F p := (Fq, F\) m for \i G (0, 1). Note that the almost reiteration property [2, V.Thm.1.5.3] ensures 



E 9 



F fl+ <-+ E- 



1+6-l/q 



^ F ff - , < o~ < 6 < 6+ < 1 . 



(3.15) 



Let now (rjj , Uj , Vj ) be a sequence of solutions to dl.7t - dl.13t inRx W 9 x W g with \rjj | + 1 1 (uj , Vj ) ||xi < B, j € N, 
for some B > 0. Writing 



d a Vj - A D Vj = -fcvj + PzVjUj —. fj , Vj(0) = £Vj =: v° 



(3.16) 



it follows from the continuity of pointwise multiplication W g x W q —> W g (owing to q > n + 2 and Sobolev's 
embedding) and ( 13.1b that 

II/jIIc(j,f„) <c(5) , jeN, (3.17) 
while (11.71 ). ( 12. j} , and the embedding E\ Fi/ q _ e with e > sufficiently small entail 

K||f 1/ ,_.<c(S), jeN, (3.18) 
for some constant c(B) > 0. Thus, from ( ETBl ), d3~T6b , dlTTTb . and (l318t . for j 6 N, 







lie 
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< <B) . 



^^WdFo.F^Wf^WFodada 



Therefore, from ( 12.11 ) we conclude that (vj) is bounded in Fi_ E . Since (/j) is bounded in C(J, Fo), we deduce from 
d316t and H II.Thm.5.3.1] that is bounded in C e ( J, i*i_2 £ ) for some e > sufficiently small. Now, taking [25, 
Thm.5.3.4,Thm.5.4. 1] into account which guarantee 

E^ 1/q = (Z?(A C )^(A 2 D )) 1 _ 1/ ^ ^ Wft 9 ^ , 

with D(A^j) denoting the domain of the fc-th power of Ad equipped with its graph norm, we obtain 

for e > sufficiently small by Sobolev's embedding theorem since q > n + 2. Consequently, 

(«,■) isboundedinC e (J,C 2+£ (n)) n W 9 . (3.19) 

But then, since 

<9 a u.,- - A D ((1 + jVj)uj) = —ctiu? - a 2 UjVj , u 3 -(0) = rfrZTj , (3.20) 

we similarly conclude that (uj) is bounded in C £ ( J, C 2+£ (Ct)), where the analogue of (13 . 14b holds due to ( 13.191 ) 
and J2, II. §5.1]. Finally, these observations warrant that the sequence (AnVj) is bounded in C e ( J, C e (f2)) while 
(-0iv? - l3 2 v 3 Uj) is bounded in C E (J, C 2+e (Cl)). From d3T6t we derive that (d a Vj) is bounded in C £ ( J, C e (fJ)) 
and similarly we derive this for (d a Uj). □ 

Noticing that C £ (J, C 2+£ (fl)) embeds compactly in C E (J, C 2+E (0,)) and C 1+E (J,C e (Cl)) in C 1+E (J, C E (fl)) 
for e £ (0, e), we deduce: 

Corollary 3.4. Any bounded and closed subset of {(r),u,w) £ R x W 9 x W q ; F(rj,u,w) = 0} is compact. 
Let now (77, it, w) £ R x W q x W g be fixed. We shall show that 

L := F^ w) (rj,u,w) e C(W q xW q ) 
is a Fredholm operator. To this end, we introduce, for a £ J, the operators Aij(a) £ C(W 2 D , L q ) by 



Au(a)</> 
A 12 (a)ip 
A 21 {a)4> 
A 22 (a)ip 



= — A £) ((l + 7«^(a) + ~tw(a))4>) + a 2 {v^{a) + w(a))(f> + 2a.\u{a) 

= -A D (ju(a)if) + a 2 u(a)^ , 

= -^a(«i(o) + w(a))</> , 

= —A D ip + 2/3 1 (u e (a) + u>(a))^> - /? 2 M a ) , 
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for a G J, (f), if) G and set 



Moreover, we define 



A(o) 



An (a) A 12 (a) 
A 2l {a) A 22 {a) 



(n \ h ._ (-&n((l + 7^(a))/ii) +a 2 v i {a)h 1 
W ~\ -A D h 2 + 2p 1 v i (a)h 2 



, a g J . 



h=(hi, h 2 ) G W„ 2 n x W? D , aeJ 



so that D G C(W? D xW^ D ,L q x L q ) for a G J, and we also define £[77] G £(W 9 xW„lV 2 fl x W 2 D ) by 



%]2 := 



z = (</>, 7/0 G W 9 X 



<q ■ 



It then readily follows from ( 13. 12b that, given z = (<fr, ip) and h = (hi, h 2 ) in '¥.',, 
equivalent to 



, q ^ q , the equation Lz = h is 

<9 Q z + A(a)z = d a h + B(a)h , aeJ, z(Q) = £[rj\z + h(0) . (3.21) 
In the sequel, we use the notation 

X := L,xL,, Xi := W,xW,, Xg := W 2e D x W 2 q % , 6 G [0, 1] . 

Let us first observe that 

Remark 3.5. The space Xi can fee equipped with an equivalent norm, which is continuously differentiable at all 
points except zero. 



Proof. According to 12211 . since Xi 



q x w, is separable, the statement is equivalent to say that the dual space 

/ 1^ ilon^i" 1 in i. \hr- ^oniiniM^ ^n:ipo IT / 



\ = W q x of Xi is separable. But, since W g is dense in L 9 , the separable space U q = L 9 / is dense in W' 9 , 



q q 

where l/q + l/q = 1, So X' x is separable. 

Investigation of ( 13.211 ) requires the following information on the involved operators: 



□ 



Lemma 3.6. The above defined operators (d a + A, 70) and (d a + D, 70) both belong to Cis(*Ki , Xn x Xi_i/ q ), and 
£[r)] belongs to /C(Xi, Xi-i/q). 

Proof. Writing 

Ai 2 (a)ip = -AD(ju(a)if) + a 2 u(a)tp 

= -div x (lu(a)V x ip) + {a 2 u(a)ip - div x (ip~/V x u(a)) } 

and using ( 13. U . it is readily seen that A can be written in the form 



Ai + A 2 



"All 


A12 




' 


A 12 





A 22 _ 


+ 


A 2 i 






with 



Recalling 



M G C(J,£(Xi,X )) , A 2 G i,(J,£(Xi_i /g ,X )) 
1 + j(v^(a, x) + w(a, x)) > 1/2 , (a,i) £ Jxfi, 



(3.22) 



due to the positivity of and u> G W 9 , it follows as in (13.6b that An(ao) and ^22(00) have maximal ^-regularity 
for each fixed ao G J. Consequently, the problem 

9 a 2i + An (a )2i +Ai 2 (a )^2 = /1 (a) , aGJ, Zi(Q) = z" , 
d a z 2 + A 22 (a )z 2 = / 2 (a) , aGJ, z 2 (Q) = 7 



admits for each / = (/i, /a) G Xo and z° — (zj, £ Xt_i/g a unique solution z = (z\, Z2) £ Xi given by 

Zi = (9 a + An(a ),7o) 1 (/i - A12 (a )z 2 ,z%) , 

Z2 = (9 a + ^22(ao),7o) {h,A) > 

and there is some constant c independent of / and z° such that 

iNik^cdi/ik + n^iix^j. 

Therefore, (d a +Ai(ao), 70) £ £is(Xi,XoxX 1 _ 1 / g ) for each ao € J, whence (<9 a +A, 70) G £is(Xi, X xX 1 _ 1 / q ) 
by (13.221 1 and Lemma [3~T1 Analogously we deduce the statement on (d a + B, 70). Since W q D embeds compactly in 

p , the assertion on £[rj\ £ £(Xi,Xl) is immediate. □ 

Based on Lemma [331 we have 

E := (<9a + A, 7o)~ 1 £ £(X x Xi) and Q := [w ^ i[rj\ (S(0, to))] G /C(X 1 _ 1 / B ) . 

We now show that L is indeed a Fredholm operator. The proof is along the lines of lEijl Lem.2.1]. 

Proposition 3.7. Lef (77, u, w) G K x W 9 x W g ant/ L = F^ u ^ w ){rj, u, w) G £(Xi). T/ien i ;s a Fredholm operator 
of index zero. More precisely, 

rg(L) = j/i G Xi ; fc(0) + %](*W + Bh, 0)) G rg(l - Q )} (3.23) 

« closed in Xi anaf 

ker(L) = {£(0 G ker(l - Q )} 

with 

dim(ker(L)) = codim(rg(L)) = dim(ker(l — Qo)) < 00 . 
Proof. Owing to ( 13.211 ) and Lemma [3~6l for z,h G Xi, the equation Lz = his equivalent to 

z = T,(d a h + Bh,0) + E(0,z(0)) , (3.24) 
(1 - Qo)z(O) = ^[97] (E(d a ft + B/i, 0)) + h(0) . (3.25) 

If 1 belongs to the resolvent set of Qo £ ^C(^i-i/<j)> then (13.241 ). ( 13.251 ) entail a trivial kernel ker(i). Moreover, in 
this case, for an arbitrary h G Xi, there is a unique z(Q) G X 1 _ 1 / q solving (13.25b , thus the corresponding z G Xi 
given by ( 13.241 ) is the unique solution to Lz = h. This easily gives the assertion in this case. 

Otherwise, if 1 is an eigenvalue of Qo £ fc(Xi-i/q)t then ( 13.24b . (13.251 ) yield the characterization of ker(L) 
and rg(L) as stated. In particular, since £ is an isomorphism, we deduce dim(ker(L)) = dim(ker(l — Qo)) which 
is a finite number because 1 is an eigenvalue of the compact operator Qo. Moreover, rg(L) is closed in Xi since 
M := rg(l — Qo) is closed by the compactness of Qo and due to Lemma l3~6l and (13. U . To compute codim(rg(L)), 
note that 

codim(M) = dim(ker(l — Qo)) < 00 , 

hence M is complemented in Xx_i/ q leading to a direct sum decomposition X\_x/q = M © N. Denoting by 
Pm £ C{Xx_xi q ) a projection onto M along N, we set 

Vh:=A(d a h + m,P M h(0)-(l-P M )iM(^(dah + Bh,0))) , h £ X x , (3.26) 

where A := (d a + B, 7o) _1 G £(X x X 1 _ 1/q , Xi), and obtain V G C(Ki) from Lemma 1331 Since 

(d a + B) {Vh) = d a h + Bh , j (Vh) = P M h(0) - (1 - P M )l[rj\(E(d a h + Bh, 0)) , 

the characterization ( 13.231 ) actually implies that V maps Xi into rg(£). Furthermore, if h G rg(L), then ( 13.231 ) also 
ensures 

Vh = A(d a h + Bh, h(0)) = h , 
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so 7 , (rg(L)) = rg(L). Thus V 2 = V with rg('P) = rg(£) is a projection and Xi = rg(L) © ker('P). Since A is an 
isomorphism, we obtain 

ker('P) = {h G Xi ; d a h + O/i = 0, h(0) G A 7 } , 
from which we deduce the equality of the dimension of N and ker('P) and thus the statement. □ 

Corollary 3.8. For k G (0, 1) and (7?, it, w) G R x W 9 x W,, 

(1 - fe)F ( „ iTO )(?7,0,0) + kF {U}W) (r],u,w) G £(Xi) 
;s a Fredholm operator of index zero. 
Proof. Since, by (EHI, 



(3.27) 



the operator 

(1 - k)F( u>w ) (t], 0, 0) + kF {u<w) (i], u, to) 
has the same structure as Fr UlW ) (rj, u,w). □ 
It follows from Lemma [331 that the operator A%, given by 

At(a)<f>:=-A D ((l+V>d a ))<l>)+ a * v e( a )<l>> aeJ > 4>£Wl D , (3.28) 

belongs to C £ (J, C(W q D , L q )), while the positivity of ensures that — A^(a) is for each a G J the generator of an 
analytic semigroup on L q . Consequently, it generates a parabolic evolution operator 11^ (a, a), < a < a < a m , in 

view of J2, II.Cor.4.4.2]. Note that IL A( (a, 0) for a > is strongly positive on W^ /q , see e.g. J3, Cor. 13.6]. We 
then set 

G(.:= I 6(a) IL^ (a, 0) da (3.29) 
Jo 

and obtain from ( 12. U and the compact embedding of in W^ 279 that G 5 G /C(W^ 2/9 ) is strongly positive. 
Thus, by the Krein-Rutman theorem, r(G^) > is a simple eigenvalue of G^ with an eigenvector in the interior of 
the positive cone W^^^ 9 ' + . Let then 

Vo -=Vo(0 ■■= >0, ker(l - r l0 G 5 ) = span{$ } , $ G int(W^ 2/<? ' + ) . (3.30) 

We define a 

<^*(a) :=II^(a,0)$o , aGJ, := / 6(a)^(a) da , (3.31) 

Jo 

and, using the notation of Subsection l3.ll 



V>*(a) := n [ 2^ 1 « 5 ](a,0)*o + / Hp^] (a, a) (/3 2 i! 5 (o-)0* (a)) dcr , aGJ, (3.32) 

Jo 

where 

#o :=e(l-e^[ 2ft , 5 ]) _1 (J b(a)J^ n [2/3l , 5] (a,a)(/3 2 ^(a)^(a))dada^ . 

Note that ^0 is well-defined since 1 — £,H[2/3 lVi ] is invertible owing to Lemma [3~2l (13.5b . and G W+ which ensure 
r(£-ff[20 lt v]) < 1- Also note, from ( 13.21 i and ( 13.6b , that 0* and ^* both belong to ' 



Lemma 3.9. The kernel of F/ UtW -)(r)o, 0,0) is spanned by (</>*,?/;*), anaf i 7 ^ ( Ujto )(?7o, 0, 0)[0*, -0*] does not belong to 
the range of F^ u w ) (rjo, 0, 0). Moreover, $0 = ^o*! 5 *- 
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Proof. Observe that (cj), e Xi belonging to the kernel of Ff u<w ) (770, 0, 0) is equivalent to 

d a (/) - A D ((1 + jv £ )<j>) + a 2 v s <p = , 0(0) = r ?0 $ , 

d a i>-A D i> + 2f3 1 v^ = (3 2 v i( f>, ^(0)=f¥, 

according to ( 13.27b and the definitions of Ti and T 2 . Now, the first assertion follows from ( 13.28b - (13.32b by solving for 
4> and ip. Next, suppose F v .( u ,w)(vo, 0, O)[0*, ip*] belongs to the range of Fi UiW \{r)Q, 0, 0). Then, in view of (13. 1 3b . 
(13.271 1. and the definition of Ti, there is e W g with 

0„0 + ^ = O, 0(0)= r;o$-^, 

so 0(a) = U.A t (a, 0)0(0), a E J and whence (1 — ri Q G^)(j)(0) = — Since $ = fyo^* by definition of 0* and 
( 13.30b , we conclude 

$o e ker(l - 770Q) n rg(l - r? Gf) 
what is impossible since ryoG^ is compact with simple eigenvalue 1. □ 



3.3. Proof of Theorem \2.2\ 

Having established the necessary auxiliary results in the previous subsection, we are now in a position to prove 
Theorem l2.2l bv applying [23, Thm.4.3,Thm.4.4]. Recall that, writing (77, u, v) = (77, u, v^ + w), the solutions (77, u, v) 
to ( 1 1 .7b -( fTTT3l > are obtained as the zeros (77, u, w) of the smooth function F defined in ( 13.1 lb . Also recall that 770 = 
770 (£) is given in ( 13.30b . 

As in the second part of the proof of Lemma |3.9l 



ker(F {UtW) (770, 0, 0)) n rg(F {u<w) ( Vo , 0, 0)) = {0} , 

whence 

Xi = span{(0*,<0*)} © r g(i J («,u))('7o,0,0)) 

by |Hl Lem.2.7.9] and Lemma [3~!9l In view of Proposition ^ .7 l and Lemma [3~!9l we may apply i23l Thm.4.3]. Therefore, 
there are e > and continuous functions 

7, : (-e, e)-*R, {61,62) : (-e, e) -> rg(F (UjW) (m, 0, 0)) 

with 77(0) = 770 and (61, 2 )(0) — (0, 0) such that the solutions to (ll.7b - dl.13b near (770, 0, v^) are exactly the semi- 
trivial ones (77, 0, Uj), 77 > 0, and the ones lying on the curve V := T + U T_ U {(770, 0, Wf)}, where 

T± := { (?7(s), s0* + s0i(s), u 4 + sip* + s6 2 (s)) ; < ±s < s} . 

Moreover, V is contained in £*, which is a connected component of the closure of 

S := {(77, u, v^+w); F(r], u,w)=0, (u, w) ^ (0, 0)} . 

Being merely interested in positive solutions, we first note: 

Lemma 3.10. The curve T + lies in R+ x W+ x W+ 

Proof. Letw s := s0» + s6»i(s) and v s := v^ + sip* + s9 2 (s). Thenu s (0) = s$ + o(s) andu s (0) = £V£ + s\& + o( s ) 
in W^ /q as s -> 0+ by (13J. Thus, it follows from $ , ^ € m^W"^ 9 '*) that u s (0), v s (0) 6 mt(T^, 2/9,+ ) 
for s G (0, e) with e > small enough, whence u s ,v s 6 W+ for s e (0, e) due to the parabolic maximum principle 
flThm.l3.5]and<0-CL13). □ 

Now, invoking Corollary 13. 41 Remark [331 and Corollarv l3.8l we obtain from B23L Thm.4.4] (see also B231 Rem.4.2.1]) 
further information about the global character of the continuum. More precisely, if £ + denotes the connected compo- 
nent of \ T_ containing T+, then £ + satisfies the alternatives: 
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(i) £ + intersects with the boundary of M x W 9 x W 9 , or 



(ii) £ + is unbounded in R x W 9 x W g , or 

(iii) £+ contains a point (77, 0, v{) with 77 ^ 770, or 

(iv) £ + contains a point (77, u, tj^ + w) with (u, w) ^ (0, 0) and (u, w) G rg(F(„ iW j (770, 0, 0)). 

Due to Lemma IXTOl the continuum £ := £ + n (K+ x W+ x W+) of solutions to (TTTTb- dl . 1 3b contains the 
curve r+. Furthermore, we have: 

Lemma 3.11. £ \ {(%, 0, v £ )} C R+ x W+ x W+ is unbounded. 

Proof. We first show that £ + does not reach the boundary of M + x W+ x W+ at some point (77, u, v) ^ (770, 0, v^). 
Suppose otherwise, i.e. let there be a sequence (rij,Uj,Vj) in £ n (K + x W+ x W+) converging toward some point 
(77, u, v) 7^ (770, 0, v^) not belonging to K + x W+ x W+. Since Uj and Vj are nonnegative, the limits u and v are as 
well. So u = or v = because (77, u, v) £ K + x W+ x W+. We claim that neither is possible. Suppose first that 
both u and v identically vanish. As Vj G W+, ipj :— Vj/\\vj\\w q is well-defined in W+, has norm 1, and 

d a ipj - A D ipj = -PxVjipj + MjUj , ipj(0) = • 

The proof of Lemma [33] shows that (ipj) is bounded in C E (J, C 2+£ (n)) n C 1+E (J, C e (fj)) for some e > and so 
we may assume without loss of generality that (ipj) converges in W+ to some ip satisfying 

d a 4> - a d tp = , v(o) = . 

Thus rp(a) = e aAD ip(0), a £ J, and ip(Q) = £H [0] ip(0) implying £r(H [0] ) = 1 by the Kr ein-Rutman theorem. 
However, this contradicts ( 13.4b and £ > 1. Next, assume u vanishes identically but v ^ 0. Then (77, u, v) — (77, 0, v) 
and the uniqueness statement of Theorem l2. ll implies v = v$. Thus, (77, 0, vg) € ^2 is a bifurcation point to positive 
coexistence states. By Lemma 1331 we may assume (vj) converges to v$ in C e (J, C 2+£ (J7)) fl C 1+e (J, C £ (ft)) for 
some e > 0. Moreover, as above we may assume that (<fij), defined by <pj := Uj /\\uj\\w q , converges in W+ to some 
<p satisfying 

d a (f>- A D ((1 + 7 ^)0) = -a 2 0v S , 0(0) = 77$ . (3.33) 

Therefore, 0(a) = (a, 0)0(0), a G J, and 0(0) = r/G^0(O). Thus 7/ = 770 by the Krein-Rutman theorem and 
(I3.301 l. This yields the contradiction (r),u,v) — (770, 0,7^). Finally, suppose v = but w ^ 0. Then we have 
(77, it, u) = (77, u, 0) what gives u = u rj with 77 > 1 by Theorem l2 . 1 1 since u G W+, and so (77, u, v) = (77, u v ,0) G Q3i 
is a bifurcation point to positive coexistence states. As above we may assume that (ipj), given by ipj := Vj/\\vj \\-w q > 
converges to some ip G W+ satisfying 

3 a iP - A d i/j = Mu v , iP(0) = ^ . 

This readily implies 1 = £r(Ht_p 2U 1) what is impossible since £ > 1 and 1 = r(-ffroi) < r(H\_p 2U 1) according to 
(t3~4b and Lemma [3~2l 

Consequently, £ + intersects with the boundary of R + x W+ x W+ only at (770, 0, vg), whence £ = £ + . So neither 
alternative (i) nor (iii) above is possible. Suppose (iv) occurs, and let (0, ip) G Xi and (77, u, + w) G £+ be with 

(u,w) = F( U)W) (rio,0,0)[<p,ip] . 

Then - u = Ti (0,770$) by (13.27b . Recall, from the definition of 0* and Lemma l3~9l that 0* = 7\ (0,770$*) 
with G int(VK g 2 D 2/9,+ ). The latter implies K77 $* + 0(0) - tt(0) G mt(W^ D 2/q ' + ) for some « > 0. Defining 
p := k0» + — 7i G W g , we obtain p = Ti(0, 770 (k$* + $)), that is, 9 p + Ajp = withp(O) = 770P + 770C/. Hence 

(1 — 77oGj)73(0) = rioU. Since m G W+ by assumption and thus U G (W^ ^ q ' + ) \ {0}, this last equation does not 
admit a positive solution p(0) according to ( 13.30b and [1, Thm.3.2] in contradiction to p(0) G mi(W^^ q ' + ) by the 
choice of k. So (iv) is impossible as well, and we conclude that £\{(?7o, 0, v^)} C K + x W+ x W+ is unbounded. □ 
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To finish off the proof of Theorem l2.2l we merely have to remark that (770, 0, v^) is the only bifurcation point. 

Lemma 3.12. There is no other bifurcation point on or on Q3i to positive coexistence states. 

Proof. Exactly the same arguments as in the first step of the proof of Lemma 13.1 II show that there is neither a 
bifurcation point (77, u v ,0) £ *Bi nor (77, 0, V% ) £ 032 to positive coexistence states. □ 

4. Proof of Theorem|23] 



As the proof of Theorem 12. 3 1 is similar to the one of Theorem 12.21 we merely sketch it and point out the necessary 
modifications. 

We shall derive a bifurcation from the branch 03 1 by linearizing around a point (77, u v , 0) with a suitable 77 = 771 to 
be determined. First note that the smooth branch!! := {(77, u v ); 77 > 1} in (1, 00) x of solutions to ( 12.3b provided 
by Theorem l2. ll extends to a smooth branch 1L := {(77, u n ); 77 > 77*} in (77*, 00) x W q passing through (77, u) = (1,0), 



where 77* £ (0, 1) and —u v £ WJ for 77 £ (77,, 1). Indeed, application of [26, Thm.2.4], 1271 Prop. 2. 5] (with e < in 



II271 Eq.(2.17)], see the proof of OOl Prop. 3.4]) shows that the branch il of positive solutions extends smoothly with 



a branch {(77(e), u^( e )); — £0 < £ < 0}, where — it^jvj £ W+ for e £ (—£0, 0). Thus, fixing e £ (—£0, 0), it follows 



that w := — u^( E ) £ W+ satisfies 



d a w — Adw = ol\w 2 , tvj(0) = rj(e)W , 



whence w(0) = r](e)H[_ aiW ]w(0) and thus rj(e)r(Ht_ aiW ]) = 1 by the Krein-Rutman theorem. Due to Lemma [3~l2l 
and (13.4b . we have r(Hi_ aiW i ) > r(H\ot ) = 1 and so 77(E) < 1. We thus get the desired smooth extension H* of il by 
choosing 77, sufficiently close to 1. Consequently, the solutions (77, u, v) = (77, u v — w, v) to problem dl.7b - dl.13b can 
be obtained as the zeros (77, w, v) of the smooth map T : (77*, 00) x W q x W q W q x W 9 , defined by 



(4.1) 



n _ f w — T( — Ao(^v(urj — w)) — 2a\Ur(w + a\w 2 + a2(u v — w)v , t]W) 

where the set W q is as in Section|3]and 

T := (d a - Ajj^o)- 1 £ C(h q x W^ D Vq ,W q ) . 

Clearly, JF(i], 0, 0) = for 77 £ (77*, 00) and the Frechet derivatives at (77, w, v) are given by 

<q - T( - A D (^(u n - w)) + A D {^v4>) - 2a\{u n - w)<j>^ 
^(w.v)^, w, v)[(f>, t/;] = I +a 2 ip(u r , - w) - a 2 vcj), 77$) | (4.2) 



and 



T , ,(77 u w)\d> M (-n-^H]- 2a <<t> + a ^ u n>*)\ 



(4.3) 



for ((f), ip) £ W 9 x W g with dashes referring to derivatives with respect to 77. It is then straightforward to modify the 
proofs of Lemma [3~!6l and Proposition l3.7l in order to derive the analogue of Corollarv l3.8l 

Lemma 4.1. For k £ (0, 1) and (77, w, v) £ (77*, 00) x W g x W q , 

(1 - k) T{ w , v ) (77, 0, 0) + k^ WjV )(r], w, v) £ £(Xi) 
is a Fredholm operator of index zero. 
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To determine the bifurcation point, let us observe that r{H^_p 2U ]) > 1 is a strictly increasing function of 77 on 
(1, 00) according to Theorem 12. II and Lemma [3~2l Since u v depends continuously on 77 in the topology of W q by 
Theorem l2.ll we obtain from [2, II.Lem.5.1.4] that the evolution operator Ut_p 2U 1 (a, 0) and hence H[_p 2U -1 depend 
continuously on 77 with respect to the corresponding operator topologies. Together with the fact that the spectral radius 
considered as a function /C(W^ 2 ^ 9 ) — > R + is continuous (see J3, Thm.2.1]), we conclude that 

(77 !-)■ r(H[_p 2Uv] )) G C((l, 00), (1, 00)) is strictly increasing (4.4) 
with lim^i r(-ff[_ i g 2Mr) ]) = 1. Defining S G [0, 1) by 

1 



lim r(H[_p 2U ]) ' 

it follows that for any £ G (<5, 1) fixed we find a unique 771 := 771 (£) > 1 with 

1 



(4.5) 



(4.6) 



We may then choose *i G int(W 2 J 2/ ' 9 ' + ) spanning ker(l - £H[-p 2Uvi ]). Define G W g x W+ by 

V>* :=n h( g 2Ui)i] (-,0)*i and 

<A* := n [2QlUiji] (-,0)$i +7V^* , $1 :=ryi(l -77i77 [2qiU7ji] ) _1 ^ 6(a)(i\ty*)(a) do^ , 

with a 

(JV^*)(a) := / n [2ai „ ](a,<r) ( - A £ ,(7u m (<j)V'*(o-)) + a 2Uj 7 i(o")V'*(o")) do- , aGJ, 
Jo 

where the invertibility of 1 — rjiH\^ aiUrii ] is due to (13.5b . The analogue of Lemma |3~9l then reads: 

Lemma 4.2. The kernel of J- \ w . v )(jji, 0,0) is spanned by (0*,^*) ant/ F^cw^} {Vii 0, 0) [0*> "0*] does not belong to 
the range of Ti WiV ^ (r]i, 0, 0). 

Proof. That ker(J r ( UI lt ,)(j7i, 0, 0)) = span{(</>*, V 1 *)} follows as in the proof of Lemma |3~9l To check the transver- 
sality condition, suppose J-^Jw^iVli 0, O)[0*, ^?*] belongs to the range of F( WtV ){T)x, 0, 0). Recall (14.211 . (14.31 l and let 
w G W q be such that u - T(f3 2 u m v,£V) = -T{(3 2 u' m 1/)*, 0). Choose r > with r*i - v{0) G intlW 2 ^ 279 ^). 
Since ip* — T(f32U Vl ip+, ^^ , *), it follows that p := rip* — v satisfies 

d a p - A D p - f3 2 u m p = PW^* , p(0) = C p , 

from which we deduce 



(l-£H [ _ f}2u ] )p(0) = £/3 2 b(a) tt hfj2Uvi] (a,o-)(u m (a)^(a)) dado, . 

Jo Jo 

However, invoking [[H Thm.3.2] and ( 14. 6h this is impossible since p(0) G int(VF 2 ^ q,+ ) by the choice of r and since 
the right hand side is positive and nonzero due to ( 12. U and the positivity of t/j* and of u' stated in Theorem l2.ll □ 



As Corollary 13.41 holds also for T, we may proceed as in Subsection [33] to derive from II23L Thm.4.3,Thm.4.4] 
that a continuum & + in (77*, 00) x W q x W 9 of solutions to d 1 .7b - (l 1.13b bifurcates from (771, u, n , 0) satisfying the 
alternatives: 

(i) 6 + intersects with the boundary of (77,, 00) x W q x W q , or 

(ii) & + is unbounded in (77,, 00) x W q x W q , or 
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(iii) & + contains a point (77, u v , 0) with 77 ^ rj±, or 

(iv) 6 + contains a point (77, u v — w, v) with (w, v) ^ (0, 0) and (w, v) G rg(Jv TO)0 ) (771 , 0, 0)). 
Moreover, near the bifurcation point, (5 + is a continuous curve 

T + = {(?7(s),m^( s ) - - s0i(s)iSip+ + s9 2 {s)) ; < s < e} 

for a continuous real-valued function r\ and some continuous W 9 -valued functions 9j with 77(0) = and 0j(O) = 0. 
Since u^( s )(0) and t/>* = ^1 both belong to uA{W^~j^^ q ' + ), it follows as in Lemma [3. 101 that T + is a subset of 
(1, 00) x W+ x W+ for e > sufficiently small. For the continuum (5, given by 

6 := 6+ n ([77*, 00) x W+ x W+) , 

we have: 

Lemma 4.3. 6 \ {(771, M^O)} is a subset of [77*, 00) x x consisting of coexistence states (ji,u,v) to 
( ll.7t - dl.13t . The continuum & is unbounded or it connects u Vl , 0) to a solution (77*, u, v) of dl.7t - dl.13t with 

u,v G W+. 

Proof. First suppose & + \ {(771, u Vl , 0)} does not reach the boundary of (77*, 00) x W+ x W+, so & — & + . Then 
neither (i) nor (iii) above is possible. Suppose (iv) occurs. Then there are (77, u v —w,v) G & and (0, ip) £ W q x W 9 
such that (w, v) = J-( w . v )(vij0, 0)[<fi, if>]. Hence p :— Kip* + ip — v G W q , with k > chosen such that p(0) belongs 

to G int(W 9 2 ^ 2/Q, ' + ), satisfies 

<9 Q p - A D p - /3 2 u m p = l3 2 u m v , p(0) = £P + , 

so that 

(l-efr [ _ ftttfl ] )p(0) = ey + f/3 2 / 6(a) / ^l-i3 2 u vl ](a,o-)(u m (<i)v(a)) dada . 

Jo Jo 

Since tj G W+ by assumption, this last equation does not admit a positive solution p(0) according to 111 Thm.3.2] in 
view of d4.6t . However, this contradicts p(0) G mt(W q i ~ r f^ q ' + ). So (iv) is impossible as well, and we conclude that 
if 6 + \ {(771, 0, v^)} does not reach the boundary of (77*, 00) x W+ x W+, then © = & + is unbounded. Otherwise, 
suppose 6 \ {(?7i, 0, v%)} reaches the boundary of (77,, 00) x W+ x W+ at a point (77, u, v) ^ (771, u ni , 0) and choose 
a sequence (r/j, Uj,Vj) in 6 (~1 ((77*, 00) x W+ x W+) converging toward (77, u, v). Since Uj and Vj are nonnegative 
and rjj > 77,, the limits u and v are nonnegative as well and 77 > 77*. So ?i = or w = or 77 = 77*. We claim that 
necessarily 77 = 77* and u, v G W+. We proceed as in Lemma l3.HI If both u = and v = 0, then the limit ?/> G W+ 
of i/jj := f/Hfj ||wg satisfies d a i/j — A^t/j = with ip(0) = £$ leading to the contradiction 1 = £r(_ff[ ]) = 6 due 
to ( 13.4b . If 77 = but 7i ^ 0, then w G W+ satisfies d a u — Arju = —a\u 2 and u(0) — r]U and thus u = u n with 
necessarily 7/ > 1 by the uniqueness statement of Theorem l2.ll Hence ip G W+ satisfies d a ip — Aijip = (32U v ip with 
ip(0) = giving the contradiction 77 = 771 by ( 14.6b . If u = but w ^ 0, then 7; G satisfies 9 a w — A^w = — /?it; 2 
and 7j(0) = ^ what is impossible according to Theorem l2. 1 I since £ < 1. Therefore, neither u = nor 77 = and we 
conclude 77 = 77*. This proves the claim. □ 

Lemma 4.4. There is no other bifurcation point to positive coexistence states on Q3i. 

Proof. The assumption (77, u v , 0) G *Bi being a bifurcation point to positive coexistence states corresponds to the 
case 77 > 1, 7i ^ 0, and v = in the proof of Lemma l4.3l and analogously implies 77 = 771. □ 

This completes the proof of Theorem l2.3l 
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5. Proof of TheoremE^l 



Again, the main part of the proof of Theorem 12.41 is a straightforward modification of Section [3] and we thus 
omit details. Let 77 > 1 be fixed. Linearization around (£, u v , 0) S T2 entails the existence of a continuum IH + in 
K x W 9 x W q of solutions to d 1 .7b - (l 1.13b bifurcating from (£ , u v , 0), where £ :— £o(v) S (0, 1) is given by 

£0 := T~lj \ ■ (5-D 

r ( fl [-fc«,]) 

Near the bifurcation point (£0, u v , 0), 1H + is a continuous curve in M. + x x and it can be shown exactly as in 
Lemma l3~TT1 or Lemma|43]that <K := n (M+ x W+ x W+) satisfies the alternatives: 

(a) <K \ {(£_ , u v , 0)} is unbounded in E+ x W+ x W+ or 

(b) D\ reaches the boundary of K + x W+ x W+ at some point (£, u, v) ^ (£ , u n , 0) with u = or v = 0. 

If (b) occurs, choose a sequence (£j , ztj , vj ) in $H n (M + x W+ x W+) converging toward (£,it,v) 7^ (£,o,u v ,0). 
Putting 4>j := Uj/\\uj ;||w 5 and ^ := Uj/Huj |]w„> we ma y assume that <j>j — »• qb and ipj — > ip in W 9 . If both u = 
and v = 0, then S W+ satisfies <9 a </> — Au0 = with </>(0) = 77$ and so 1 = ?/r(i?[ ]) contradicting ( 13.41 ) and 
?7 > 1. If u ^ but v = 0, then u — u v according to Theorem l2.ll and ip G W+ satisfies 9 a ?/> — A^^ = f32U v ip with 
■0(0) = ^W. Hence £ = £0 by J5.lt what is impossible since (£, u, w) 7^ (£0, u v , 0). Therefore, the only possibility 
is u = but v ^ 0. In this case necessarily £ > 1 and u = in view of Theorem l2.ll So 91 joins up with Ti at 
(£i 0) We remark that then the relation 

V r(Gs) = 1 (5.2) 

with given in (13.29b must hold, since qb 6 W+ satisfies ( 13.33b . That no other bifurcation point(s) on T2 or Ti 
exist(s) is immediate by the previous observations. This yields Theorem l2.4l 



Remark 5.1. Since the operator in (13 .28b does not yield a suitable maximum principle, there is no analogue of 
Lemma [Ol f or the spectral radius of and the only information we have on r{G^) is that it is positive for each 
£ > 1 as observed in Section 13.21 Consequently, given r\ > 0, we cannot decide a priori whether (15.2b holds for 
some £ > 1. However, for i] > 1, if VI joins up with %\, then ( 15.2b must occur and the connection point (£, 0, v^) on 
Ti is determined by this relation. Then again, the relation (15.2b is also a sufficient condition for the existence of a 
continuous curve of positive coexistence solutions bifurcating from %\. 

The same difficulty arises when considering bifurcation from Ti with respect to £ when rj < 1 is fixed. In this case, 
(15.2b is again a necessary and sufficient condition for the existence of a bifurcation point on Ti to a curve of positive 
coexistence states, which then extends to an unbounded continuum. 
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